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Abstract 


Wave kinetic theory has been developed to describe the statistical dynamics of weakly nonlinear, dispersive waves. 


However, we show that systems which are generally dispersive can have resonant sets of wave modes with identical 
group velocities, leading to a local breakdown of dispersivity. This shows up as a geometric singularity of the resonant 
manifold and possibly as an infinite phase measure in the collision integral. Such singularities occur widely for 
classical wave systems, including acoustical waves, Rossby waves, helical waves in rotating fluids, light waves in 
nonlinear optics and also in quantum transport, e.g. kinetics of electron-hole excitations (matter waves) in graphene. 
These singularities are the exact analogue of the critical points found by Van Hove in 1953 for phonon dispersion 
relations in crystals. The importance of these singularities in wave kinetics depends on the dimension of phase space 
D - {N - 2)d (d physical space dimension, N the number of waves in resonance) and the degree of degeneracy 6 
of the critical points. Following Van Hove, we show that non-degenerate singularities lead to finite phase measures 
for D > 2 but produce divergences when D < 2 and possible breakdown of wave kinetics if the collision integral 
itself becomes too large (or even infinite). Similar divergences and possible breakdown can occur for degenerate 
singularities, when D - d < 2, as we find for several physical examples, including electron-hole kinetics in graphene. 
When the standard kinetic equation breaks down, then one must develop a new singular wave kinetics. We discuss 
approaches from pioneering 1971 work of Newell & Aucoin on multi-scale perturbation theory for acoustic waves 
and field-theoretic methods based on exact Schwinger-Dyson integral equations for the wave dynamics. 

Keywords: wave kinetics; resonant manifold; Van Hove singularity; turbulence; quantum Boltzmann equation; 
graphene. 


1 . Introduction 

Wave kinetic equations were first introduced by Peierls in 1929 to discuss thermal transport by phonons in crystals 
m and have since been extended to a large number of quantum and classical wave systems an I a a. However, 
many questions remain about the validity of these equations, in which wavenumber regimes they hold and under 
what precise assumptions on the underlying wave dynamics. Some authors based their derivation of the wave kinetic 
equation on the RPA (“random phases and amplitudes”) assumption aa , but it has been cogently argued ana that one 
can relax the RPA assumption to dispersivity of the waves. In this paper we argue that the dispersivity requirement is 
more stringent than what has commonly been understood. In particular, a wave system that is generally dispersive can 
experience a breakdown of dispersivity locally in the V-wave phase space which renders the wave kinetic equation 
ill-defined. Consider, for example, a general 3-wave equation for the evolution of the wave-action spectrum, n(k, f) : 



i=(-l,i2,S3) 


(«(k2)n(k3) - i2«(k)n(k3) - i3n(k)n(k2)}. 


( 1 ) 


X 


Here we use the shorthand notations k = (k, k2, k3), s - (s, S 2 , ^3) with a-b-ab + 02^2 + 03^3. The integer s labels 
the degree of degeneracy of the waves, or the number of frequencies a>s corresponding to a given wavevector k. We 
have assumed above that s = ±1 and Ws(k) = saj(k), appropriate for systems second-order in time where there are 
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two waves traveling in opposite directions. Because of the Dirac delta functions, the collision integral is restricted to 
the resonant manifold: 

“Rk = {k2 : sw(k) + S2w(k2) + S 3 w(k 3 )|. , . = o}. 

^ Ik3=-j'3(.yk+52k2) ' 

An immediate issue in making sense of the kinetic equation is the well-definedness of the measure on Rk ■ 

f d‘‘k 2 f d‘'k 3 6(s ■ a>(k))S‘'(s ■ k) = f .. x, ■ (2) 

J J ~ ~ M - Vw(k3)| 

In the last equality we have used a standard formula to express the measure in terms of the surface area S on the 
imbedded manifold in Euclidean space (e.g. see Q, theorem 6.1.5). The difficulty occurs at critical points where 

Vw(k2) = Vw(k3) 


and the denominator vanishes. Since V(u(k) is the group velocity of wavepackets, such points correspond physically 
to wavevector triads at which the dispersivity of the wave system breaks down locally and for which two distinct 
wavepackets from the triad propagate together for all times with the same velocity. The kinetic equation may be 
ill-defined for k at which the denominator produces a non-integrable singularity. 

An exactly analogous problem was studied in 1953 by Van Hove lH for the phonon density of states in a crystalline 
solid, which is the function defined by 


g(oj) 


^L‘^y f cfk 




ui,(k)) 


y r dS(k) 

(2nY ^ Jaj,(k)=oj |Vwj(k)| 


where the the sum s is over the branches of phonon frequencies and the k-integral is over a Brillouin zone of the crystal. 
Critical points of the phonon dispersion relations, where the group velocity Va;j(k) vanishes, may lead to singularities 
in the density of states. Indeed, by a beautiful application of the Morse inequalities ||9l. Van Hove showed that the 
homology groups of the Brillouin zone (topologically a (f-torus) make such critical points inevitable. He showed 
further by a local analysis of the critical points using the Morse Lemma ||9]|Tni that the resulting singularities are non- 
integrable for space dimension d -2, giving rise to a logarithmic divergence in the density of states, and are integrable 
for d - 3, producing divergences only in the derivative (cusps). As we shall see in the following, the local analysis of 
Van Hove carries over to the corresponding point singularities in the measures Q on the resonant manifolds, which 
we propose to call resonance Van Hove singularities because of their close physical and mathematical analogies with 
the singularities considered by Van Hove. 

It is important to emphasize that the global topological argument for existence of critical points given by Van Hove 
does not carry over to the wave kinetic problem, even when there is a spatial lattice and the reciprocal wavenumber 
space is a cZ-torus. The most essential difference is that the density of states function involves all level sets of the 
dispersion relation, whereas only the zero level sets are relevant for wave kinetics. Nevertheless, topological criteria 
for existence of singularities can be sometimes exploited in wave kinetics. For example, in the absence of singular 
points, the resonant manifold can deform continuously from Rk^ to Rk^ for any ki and k 2 . It is thus sufficient to 
establish the existence of singularity if one can find certain topological invariants of two manifolds Rk^ and Rkj for 
kj + k 2 , e.g., the fundamental group or the homology group, which are different. We shall not develop general criteria 
along this line but instead show by many concrete examples below that critical points often appear on the resonant 
manifold in practice. Figure 1 plots a simple example in which physical space is a cubic lattice and wavenumber 
space is a 3-torus, which can have relevance for numerical simulations of wave turbulence on a rectangular grid (for 
more discussion, see section 2.\.2\ . This example illustrates the general geometric feature of resonant Van Hove 
singularities that the “resonant manifold” is no longer a true manifold, locally diffeomorphic to Euclidean space. 
When zero is a regular value of the resonance condition, then the Regular Value Theorem/Submersion Theorem IfTTl . 
guarantees that the “resonant manifold” is indeed a smooth manifold, but this is not usually the case when zero is 
a critical value. If the dispersion relation is not a Morse function (i.e. a smooth function with only non-degenerate 
critical points), then the singularities can occur also along critical lines and surfaces within the resonant manifold. 


* All of the resonant manifolds exhibited in this paper are plotted with MATLAB, using contour for D = 2 and Isosurface for D >3 
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Figure 1: 

1 

l+log 2 cos(l/4) 


Three-wave resonant manifold = {p : ti;(p) + <i;(k - p) - = 0| for dispersion relation <x)(k) = [42f_| sin^ (y)] 

on the 3-torus [-n, ;r]^, and for the specific wavenumber k = (0,0, 1). There is a critical point (black dot) at p = k/2. 
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The critical set on the singular manifold may have effects on the wave kinetic theory either little or drastic, 
depending on their character. As in the analysis of Van Hove, isolated critical points have only very mild consequences 
for 3-wave resonances in space dimensions d > 3, but lead to divergences for d - 1,2, which may alter the naive 
predictions of the wave kinetic theory or vitiate the theory entirely. For higher-order wave resonances and in higher 
space dimensions, critical lines and surfaces can have similar effects. We shall give below examples of resonance 
Van Hove singularities for many concrete systems, including acoustic waves in compressible fluids, helical waves in 
rotating incompressible fluids, capillary-gravity waves on free fluid surfaces, one-dimensional optical waves, electron- 
hole matter waves in graphene, etc. A pioneering study of acoustic wave turbulence by Newell & Aucoin ca argued 
that the breakdown of dispersivity of sound waves on the resonant manifold leads to a different long-time asymptotics 
and a modified kinetic equation. The breakdown for acoustic waves is extremely severe, in that all points on the 
resonant manifold are critical (for more discussion, see section 2.1.1 1 . As we shall discuss in this work, even a single 
critical point on the resonant manifold can lead to breakdown of standard wave kinetics, which then must be replaced 
with a new singular kinetic equation with a collision integral resulting only from resonant interactions in the critical 
set. 


The detailed contents of this work are as follows. In sectionj^we present many physical examples of wave systems 
with various degrees of divergence in the phase measure and collision integral due to local breakdown of dispersivity 
and we discuss some of the specific mechanisms which can produce such critical points. Next we present in section 
l^a few general results on the effects of these singularities on the local finiteness of the phase measure, including the 
argument of Van Hove for non-degenerate critical points and some analysis of the effects of degeneracy. Finally, in 
section we consider briefly the singular wave kinetics that becomes relevant when the standard kinetic equation 
breaks down, developing the multi-scale asymptotics suggested in IfT^ and also the field-theoretic approach to gen¬ 
eralized kinetic equations suggested in ns. The main example we consider in this section, of some physical interest 
by itself, is the kinetics of the electron-hole plasma in defectless graphene IHIISl. Full details of the derivation of 
the quantum kinetic equation, which are relatively standard in the wave turbulence community, are relegated to an 
Appendix. 


2. A Case Study of Resonance Van Hove Singularities 

In this section, we want to consider singularities that occur for various dispersion relations commonly discussed 
in the literature. For simplicity, we assume throughout that the dispersion relation depends only upon wavenumber, 
not position. We begin with the simplest case of 3-wave resonances. 
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2.1. Triplet Resonances 

The conditions of resonance for triplet interaction can be written, without loss of generality, as 


S2w(p) + S3w(q) = w(k) 

(3) 

S2P + S3q = k 

by taking s' = -s,s —> s,, and the condition for singularity as 

(4) 

Vw(p) = Vw(q). 

(5) 


2.1.1. Isotropic Power Law 

In wave turbulence literature, the most commonly considered dispersion relation is an isotropic power-law; 

w(k) = Cr. (6) 

Resonance requires a ^ 1 a. The singularity condition Q becomes 


and thus 




( 8 ) 

We treat first the case a > I, where clearly ([^ implies p - q and Q then implies p = q. Substituting back into 
0,0 gives 

= (S2 + S3)p'^, k = (^2 + S3)q. 

If p = q - 0 , then also k - 0, which is the trivial resonance with all members of the triad zero. If p - q > 0, then we 
can have either S 2 + = 2 or S 2 + ■*3 = 0. However, S 2 + ■*3 = 2 implies both = 2p“ and k - 2p, which requires 

a = 1 , a contradiction. If instead S 2 + ■*3 = 0, then k = 0 and any p = q are allowed. The “resonant manifold” is 
here all of space. However, for any dynamics which preserves the space average of the wave field, the interaction 
leaves the zero-wavenumber Fourier amplitude invariant. The dynamics of the k = 0 mode is then null and there is no 
interest in the kinetic equation for that wavenumber. (This conclusion is, however, a bit too glib, as we shall discuss 


further in Section 2.1.3 below.) Thus, no singularities of a nontrivial type are allowed for a > 1. 

The case a = 1 is different, as has been long understood iflhl [T^ . For this situation, the condition 0 gives no 
restriction on p, q and 0 is the condition of collinearity: 


p = q. 


(9) 


For S 2 = *3 = 1, the conditions 
are satisfied by 

For S 2 - -Si = 1, the conditions 
are satisfied by 


k - p + q, k = p-i-q 
p=,6k, q = k-p = (l-yS)k, ;ee[0,l] 
p = k + q, p = kH-q 
p = yk, q = p - k = (y - l)k, ye[l,oo) 


Finally, for -S2 = S3 = 1, we just reverse p «-> q in the last two equations. Because of our initial choice s = -1, all of 
the above solutions have p-k > 0. Taking k —> -k, |k| = k ^ | - k| = k gives solutions with p-k < 0. Note that the 
resonant “manifolds” where 


= S2w(p) + S3w(s3(k - S2P)) - w(k) = 0 


( 10 ) 
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are now straight line segments or rays, and not manifolds at all. The entire “manifold” is in the critical set where 
Vp£’^^''^(p; k) = 0. It is a consequence of the Morse Lemma that non-degenerate critical points must be isolated, but 
the critical points here form a continuum. In fact, these critical points are all one-dimensionally degenerate, since the 
Hessian V ® k) has k as an eigenvector with eigenvalue 0. 

The dispersion law (j^ with a = 1 occurs physically, for example for sound waves, and all of the above facts 
have been noted in previous discussions of acoustic wave turbulence insiiii. This was termed a “semidispersive” 
wave system, since sound waves are nondispersive on the resonant manifolds, along the direction of propagation, but 
there is lateral dispersion due to angular separation of wave packets. The usual wave kinetic equation is no longer 
applicable, as the standard phase measure on the resonant manifold becomes ill-defined. Various proposals have been 
made to derive related kinetic equations by generalizing the arguments for strictly dispersive waves lfT2l or by taking 
into account nonlinear broadening of the resonance lfT3l . The breakdown of dispersivity is quite severe for sound 
waves, but, as we have seen, it is an isolated phenomenon for wave systems with power dispersion laws, occurring 
only for a - 1. This may have led to an expectation that dispersivity breakdown and the associated singularity of the 
resonant manifold is otherwise absent. However, we shall now show by further examples that it occurs more widely, 
although in generally less severe forms than for sound waves. 


2.7.2. Lattice Regularization of Power Laws 

The work of Van Hove suggests that the presence of a spatial lattice could facilitate the appearance of resonance 
singularities. We show that this expectation is correct, by considering a “lattice regularization” of the power-law 
dispersion of the previous section. If physical space is a cubic lattice al/, then Fourier space is the li-torus = 

f ] • The dispersion law must be a periodic function on A* = [-f, - A natural replacement for is therefore 
4 Yfi=\ sin^(^), the discrete Fourier transform of the lattice Laplacian -A^ defined by 2nd-order differences. We thus 
consider 


(y(k) = C 




( 11 ) 


This dispersion law is a toy model that does not describe any physical system but that we use as the simplest example 
of the effect of a lattice. It could appear in a computational algorithm to solve a wave equation with a power-law 
dispersion relation on a regular space grid. A lattice formulation of the equations can also be useful for theoretical 
purposes as an “ultraviolet regularization” to avoid high-wavenumber divergences, a mathematically cleaner alterna¬ 
tive to the Fourier-Galerkin truncation of the wave dynamics employed in 0. Of course, lattices exist physically in 
nature, e.g. in crystalline solids, and similar singularities must exist in certain cases on the resonant manifolds in the 
quantum Boltzmann equations which describe particle transport in crystals. Since this a 4-wave resonance, we discuss 
it in more detail in subsection 2.2 below. It would be interesting to investigate whether such resonance singularities 
exist more widely for quantum transport in crystals. 

We now turn to our simple example. Without loss of generality, we take lattice constant a - 1. We first show that, 
as for the continuous case, nontrivial resonance requires a > 1 : 

Proposition 2.1. If a ^ 1, the resonant manifolds - {p : k) = 0) contain only two points 0 and k. 


Proof We introduce the notation 


|k|A. = 


d 

4 ^ sin2(k,/2) 
1=1 


1/2 


k e 27rT‘' 


for the “vector norm” on the torus. We see that this function satisfies the triangle inequality, |p + qIa* < IpIa* + IqIa*, 
with equality only for p = 0 or q = 0 (mod A*), because of an observation of lb) that |k|A» = |z(k)| for z(k) 6 


defined by z/k) = 1 - e'^‘ with z^(p -H q) 


Zy(p) -I- e‘^t^j(q). In particular note that e' 


ipj 


z/p) = A ■ Zj(q) for some real 


A 6 [0, oo] if and only if either pj = 0 or qj - 0 (mod 27r). It follows from this fact that the standard proof of lack of 
resonances for k in Euclidean space 0 carries over. For completeness we give the details. 

Consider first S 2 = ■*3 = 1, so that k = p -H q. Then, 


|k|A- = Ip + qU- < IpIa- + IqU- < [IpIa- + IqlX-] 


I fa 


for a ^ 1. 









so that w(k) = |k|", < |p|", + |q|^, = i2w(p) + with equality only for p = 0 or p = k. 

Now consider S 2 + sj, - 0, for example S 2 = = 1- In that case p = k + q and the previous argument shows that 

|k|“. +|q|X. >|plX.,or 

m(k) = |k|X. > IpIX. - IqIa. = ^2w(p) + i3t^(q), 

with inequality only for p = k or p = q. The second possibility requires k = 0 , however, and is not of interest. 
Repeating the argument for 53 = -i 2 = 1 leads to the same conclusion with 2, p <-4 3, q so that equality requires q = k 
or q = p. The first possibility is equivalent to p = 0 and the second to k = 0 , again not of interest. □ 

We now establish results on existence of critical points, separately for the cases i 2 • ■S 3 = 1 and S 2 ■ sj, - -I '■ 

Proposition 2 . 2 . For any a > I, there always exists iifor which the resonant manifold contains nontrivial critical 
points. 

Proof The singularity condition 

IpIaT^ sin(p,) = Iql^:^ sin(q',) for i = 1, • • • , t/. (12) 

is solved by either p = q = |k or by p = q = |k + ?r, where n - {n,n, Observe that the notation “k/2” is 
actually ambiguous, because changing the components of kby integer multiples of In yield different points depending 
on whether the integer is odd or even. We remove that ambiguity by requiring that -n < kj < ji, j — 1,..., c/ so that 
both jk and jk + ;r are well-defined (mod A*). In either case, p H- q = k (mod A*). 

The resonance condition a;(p) + w(q) = a;(k) then yields for q = p = jk 



where we define the average (/(k)) = Yfi=\ Note (cos^(|)) = 5 when k = ;r and and 

approaches 1 when k —» 0 . Thus for any 1 < a ^ 2 , there always exists k that solves the condition ( [T3] l by the 
intermediate value theorem. 

On the other hand, for p = q = ^k. + n, the resonance condition w(p) + m(q) = cu(k) yields 



where we now define {/(k)> = Yj‘i=i /(^i) cos^(^;/4)/[2f=i cos^(^,74)]. Note {sin^(k/4)) = 0 for k = 0 and equals j 
for k = ;r. Thus for any a ^ 2, there always exists k that solves the condition ( [l4| ). 

Note that in the borderline case a = 2, the singularity condition becomes simply 

sin(/ 7 ;) = sin(^; - p,) for i - I,-- - ,d. (15) 

and one can easily check that there are only resonant solutions for k = ;r. For this choice, every point of A* is in the 
resonant manifold, since 


• 2/P 
sm - 

'A • 2 /? 
- - 1 - sm - 

'l = sin^(^ 

h\ 2 IP 

- + cos - 

1 = sin2(^) 


! / \2 


l ) \2 

\ 2 ) 


Every point of the torus is a degenerate critical point in this case. □ 
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The condition ( |l3|l provides the example in the Introduction. For k = (0,0,1), one gets (cos^(k/4)) = cos^(l/4), 
which reduces ( |13| to — cos^(l/4), or a = 1/(1 + log 2 cos(l/4)) = 1.0477. One can likewise use condition 

(14 1 to construct examples for a > 2. For instance, if one chooses k = (3,3,3), then (sin^(k/4)) = sin^(3/4), and one 
obtains from (14 1 that a = 1/(1 + log 2 sin(3/4)) = 2.2367. This example is plotted in Fig. 2. 



Py 


Figure 2: Resonant manifold = IP : m(p) + tuCk - p) - tu(k) = 0| for dispersion relation (i)(k) = [4 (I-)] ^ ,«= ™ 

the 3-torus [jTZ, |;r] X [- jtr, jtr] X [jTT, ^k], and for the specific wavenumber k = (3,3,3). There is a critical point (black dot) at p = k/2 - n. 

We finally discuss briefly the case S 2 ■ = -1, or, without loss of generality, S 2 = 1, = -1- We present here 

no construction of critical points for all values of a > 1, but we can establish their existence in specific instances. A 
simple example is provided by exploiting a general criterion for the existence of critical points, namely, that there be 
a pair of distinct wavevectors ki and k 2 such that the topological invariants of "Rki and are different. For example, 
in the case a = 5/2, critical points must be present because there is a change of topological invariants of the resonant 
manifold as k is varied, see Fig. 3. The manifold for ki = (.2, -.2) is connected and belongs to the trivial 1st homology 
class on the 2-torus, while the manifold for k 2 = (.18, -.2) has two connected components in the same non-trivial 1st 
homology class (both winding around the same hole of the 2-torus). Thus, as the wavevector is changed continuously 
from k] to k 2 , a critical point in ‘Rk must occur at some intermediate value of k. Such a critical point is shown in the 
middle subplot of Fig. 3. 





Figure 3: Resonant manifold for dispersion relation m(k) = sin^(y)] with a = 2.5, on the 2-torus [-ff/2,3/2n'] X [ff/2,5/2n’]), 

and for the specihc wavenumber k = (.2, -.2) on the left and k = (.18, -.2) on the right. For the specific numerically approximated wavenumber 
k = (.1887713, -.2) (middle figure), there is a critical point (black dot) at p = (-1.4716,4.6016). 

The examples in this section are only toy problems meant to provide some analytical insight. In section]^ we 
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shall show that non-degenerate singularities at a single point, like those in the examples of this section, can lead in 
appropriate circumstances to an infinite phase measure on the resonant manifolds. For 3-wave resonances, the phase 
measure is finite for dimensions d > 2 but logarithmically divergent for d - 2. We see below that such non-degenerate 
singularities can arise from dispersion laws in physically relevant models. 

2.1.3. Anisotropic Dispersion Relations 

Many wave dispersion relations in physical systems are anisotropic, either power-laws or more general forms. 
Here we show by several examples that anisotropy can lead to critical points on the resonant manifolds. A simple 
class of examples are systems whose dispersion relation has the form 

«(k) = kupik) 

where k - |k| is the Euclidean norm and k\ is a component in a distinguished direction. In cases such as this, for any 
“slow mode” k with ki - 0, 


= fp : k) = p\(tp{p) - ipiq)) = 0, q' = |k - pi). 

Since 

Vp£'^^(p;k) = {ip{p) - ip{q))tx + p\[^'{p)p + 

it follows that the subset 

= {p : Pi = 0 & ^(p) = ipiq), ^ = |k - pl) c 

consists of critical points satisfying Vaj(p) = Vm(q) oc Ci. Note also that when pi = 0 the Hessian matrix becomes 

V ® V£'^^(p; k) = ip'(p)(ei ® p -r p ® ei) -h (p q) 

so that the critical points p e with (fi'ip), tf'iq) + 0 are non-degenerate for d - 2, while for d > 2 the Hessian 
matrix is degenerate. A null eigenvector of the Hessian for any d > 2 when p, q are non-parallel is given by the vector 
i/3'(p)q-‘- - where p-*- is the component of p perpendicular to q and vice versa for q-*-. When p||q, then any 

vector orthogonal to both Ci and p (or q) is a null eigenvector. For d > 3, furthermore, any vector orthogonal to all 
three vectors Ci, p, and q (or k) is a null eigenvector. 


Rossby/drift waves. A well-known example of an isotropic dispersion relation occurs for Rossby/drift waves, where 
d = 2, the 1-direction is zonal, (f(k) - -jSp^ l(\ + p^k^) with (3 the beta parameter (meridional gradient of the Coriolis 
frequency) and p the Rossby radius IfTTI . In that case, the condition y?(p) = (p{q) can be easily solved to give k-(p - 
k/2) = 0. For a general “slow” wavevector k = (0, ky), the resonant manifold is the union of the two lines Px -0 
and py - kyl2, see Fig. 4, middle panel. The point of intersection {0,kyl2) is a nondegenerate resonance Van Hove 
singularity. Note that for generic k, the manifold R^^ for the system of Rossby waves is instead diffeomotphic to a 
circle, as shown in Fig. 4, left and right panels. 

It is easy in this example to exhibit explicitly the logarithmic divergence in the phase measure dS/\VE\ on the 
resonant manifold, which is expected for d -2. A simple calculation gives 


Vp£++(p;k) 


2/3p*ky 


(1 +p^p^)(l +p^q^) 


Py-j]^ + Pxy 


on the resonant manifold for k = (0, ky), vanishing approaching the singular point. Since ^^^(p; k) oc pxipy - ky/2), 
the divergence has the general form of the integral 


J'J' dxdyd(xy) - J'J' dxdy 


n%) + 

IjcI lyl 


One might dismiss this singularity as dynamically itTelevant, since the zonal flows with kx - 0 have vanishing nonlin¬ 
earity. This may be easily verified for the Chamey-Hasegawa-Mima equation HD, 

d,(p^ AiJ/ -if/) -jSdxifr + AJ(if/,p^ Alp) -0, 
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Figure 4: Resonant manifold k = (cos 8, sin 6) for Rossby/drift waves, plotted as gray lines. The black arrow is the vector k. Here 6 = f -0.05 

for the left panel, 6 = | for the middle panel, and 0 = | + 0.05 for the right panel. 


with J{f,g) - fxgy - fygx the Jacobian, which vanishes whenever one of the functions is independent of x (or of y). 
This argument is correct, but must be made carefully. 

The delicate point is that the critical point for the resonant manifold with kx - 0 implies not only divergent 
phase measure on that manifold but also extremely large phase measures on adjacent manifolds with kx very small. 
Note for such kx that the resonant manifold locally for p near {Q,kyl2) is a hyperbola whose equation is, to leading 
order. 


(Px - 



'^)-^-kk 
2 8 


1 + 

1 + 


Near p = (0, kyl2), to leading order. 


V£++(p;k) 


--[ (3kx(2px - kx) + ky{2py - 


ky)) X + [kx{2py - ky) + ky{2p x “ ^x)) 


Therefore on the resonant manifold 


|V£++(p;k)| ^ 


(1 + 


V 


I6\kxky{ky - 2py)(kx - 2px)\ 


(1 + f ) 


6klkj 


1 + \p^kj 
1 + p^kj. 


which implies a phase measure which is 0{\kx\ ^). Fortuitously, however, all of the standard dynamical models of 
Rossby/drift waves ini have interaction coefficients vanishing proportional to \kxPxqx\'^^^ for any of the three 

wavevectors in the “slow” set. In these models, the singularity in the phase measure for small kx is cancelled by 
the interaction coefficient |//k,p,qP ~ IkxPx^xl rendering the collision integral finite even as kx tends to zero. If the 
interaction coefficient had vanished more slowly then \kx\^^^ in the limit, then the collision integral for near-zonal flows 
could become large, threatening the validity of the kinetic description imiiiiiia. These considerations apply more 
generally, e.g. to the isotropic power-law dispersion relations discussed in section : 


2 . 1.1 


Inertial waves. Another example is inertial waves, where d - 3, the 1-direction is the rotation axis, and <p(k) - 2Q./k 
with Q. the rotation rate ll20l . This case is geometrically quite similar, but extended to li = 3. For a generic wavevector 
k the resonant manifold is diffeomorphic to a sphere. For a “slow” mode with k-Q = 0, however, the resonant 


^In the case of an isotropic power-law dispersion relation for |k| <K |p| 

Vp£:+“(p; k) = q-IpF^^P - fflp - kr^^fp - k) k a\pf^^ [k*^ + (a - l)t:||], 

where k^, iy are components of k perpendicular and parallel to p, resp. A lower bound follows that |Vpii^“(p; k)| > conir.|pP“^|k| and the limit 
k —» 0 gives a finite collision integral if the interaction coefficient vanishes no slower than |k|*/^ 
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manifold is a union of two planes, one orthogonal to fl and one orthogonal to k. For example, if Q = Qz and 
k = (kx, 0,0), these are the planes Pz -0 and px = kxl2. See Fig. 5. Consistent with our general discussion for d = 3, 
the critical subset of the resonant manifold is 1-dimensionally degenerate, given here by the intersection of the two 
planes. The phase measure on the resonant manifold is logarithmically divergent in the vicinity of the singular line. 
This result shows by example that, while non-degenerate critical points produce integrable singularities for d > 2, 
line singularities can lead to divergences in three dimensions. Although geometrically quite similar to the case of 
drift waves, the situation is dynamically very different. While for drift waves the nonlinearity vanishes for the “slow” 
modes, in the case of inertial waves the “slow” modes with k-Q = 0 correspond to a strongly interacting system 
described by 2D Navier-Stokes dynamics. It has been argued convincingly that the kinetic theory for inertial waves 
must break down in the vicinity of this 2D plane of “slow” modes, as there is there no separation of time scales 
between fast linear and slow nonlinear dynamics Il20l . Here we see that there is also a breakdown in the fundamental 
assumption of dispersivity of waves. There is an infinite set of wavevector pairs of “slow” modes with identical group 
velocities along the rotation axis and triads formed from these pairs produce a diverging contribution to the phase 
measure. Thus, the kinetic equation for inertial waves is not even well-defined for wave action non-zero in the vicinity 
of the “slow” 2D modes. 



Figure 5: Resonant manifold k = (1,0,0) for inertial waves with rotation about the z-axis. The black arrow is the vector k. 


Internal gravity waves. A third example of an anisotropic dispersion law in = 3 of a slightly different sort is internal 
gravity waves, where the 1-direction is vertical (the direction of gravity), and cu(k) = Nknlk, with kn - the 

magnitude of the horizontal component and N the Brunt-Vaisala frequency ED. In this case 

Vw(k) = kx(kxei). 

knk 

For a “slow mode” with only vertical variation (kn - 0), it is straightforward to see that = {p : Ph - 0), i.e. the 
resonant manifold is the 1-axis or the set of slow modes. Since in that case p || Ci, Vcu(p) = 0 and the entire resonant 
set consists of (degenerate) critical points. As in the case of Rossby waves, however, the nonlinear interaction 
coefficient of the Euler-Boussinesq system vanishes rapidly near the set of slow modes (see ED, eqs.(62) and (63)) 
and this singular manifold is not dynamically relevant in wave kinetics. 

Another interesting phenomenon is seen in the resonant manifold of internal gravity waves when k is a 2D mode, 
with ki = 0. It is easy to show in that case that the only possible critical points in are also 2D modes and, 
because of the restriction pnlp + qnlq - T the only allowed values are p = 0 and p = k. A plot of the resonant 
manifold for k = Ci in Fig.6 below shows that geometric singularities indeed occur at p = 0, k. However, these do 
not correspond to ordinary critical points where Vp£’^^(p; k) = () but to points p instead where |Vp£’^^(p; k)| = oo. In 
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this example, the singularities are cube-root cusps, since the equation for the resonant manifold is given near p = 0 in 
cylindrical coordinates by pi = ±(2k^pH)^^^, to leading order. Although we have considered here a 2D mode k, the 
resonant manifold of internal gravity waves exhibits similar cusps at p = 0 and p = k for generic k, because of 
the divergence of Vw at the origin. 



Figure 6: Resonant manifold k = (1,0,0) for internal gravity waves with vertical direction along the z-axis. The black arrow is the vector k. 

In general, we shall use the term “pseudo-critical point” for any point p on a resonant manifold where 
^^‘^^(pik) is non-smooth in p. Although such points may give rise to geometric singularities, they do not usually 
produce an infinite phase measure. In fact, the density of the phase measure with respect to surface area (Hausdorff 
measure) vanishes at points where |Vp£’^^(p; k)| = oo, and thus the phase measure is locally hnite whenever the 
Hausdorff measure is locally finite. The latter condition may easily be checked for the pseudo-critical points in 
Fig. 6 by using the standard formula for element of surface area in cylindrical coordinates to obtain near p = 0 that 
dA = p^Jj^dpH, which has locally a finite integral. Note furthermore for the Euler-Boussinesq system that 

the pseudo-critical points are not dynamically relevant in wave kinetics, since the nonlinear interaction coefficient 
vanishes when all modes in the triad have zero vertical wavenumber ini- 

Summary. As these examples show, anisotropy—whether power-law or other type—can readily lead to critical points. 
In most of the common cases that we have examined, the singularities in the resonant manifold are protected by 
vanishing nonlinearity from having any dynamical effects. Such protection is by no means guaranteed. The case of 
inertial waves presents an opposite case, where the singularity is associated to strong nonlinearity and a breakdown 
of the wave kinetic theory. More generally, the singularities can have intermediate effects between none at all and 
complete breakdown of wave kinetics. We shall present examples of this in the next section. 


2.2. Quartet Resonances 

Resonance Van Hove singularities also occur in 4-wave systems, for which the collision integral has the standard 
form; 


Ck[n] 



d‘‘q d‘‘{ 6‘‘(p -H q - k - f)5(w(p) -h w(q) - w(k) - w(f’)) 


xn(p)n(q)«(/’)n(k) 


1 1 1 
n(0 n(k) n(p) 


1 

n(q) 


H 


k.ft 


p.qi 


(16) 


When the wave dynamics is Hamiltonian and 3-wave resonances are absent, the collision integral can be brought 
to the above form by a canonical transformation of the Hamiltonian system ||2l- The resonant manifold is now 


= {(p,q) : E^^\p,q;k) = 0), with 


q; k) = w(p) -H cu(q) - a)(k) - w(p -i- q - k). 


(17) 
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Thus can be expected to be a {2d - l)-dimensional surface embedded in a Euclidean space of dimension D - 2d. 
As a matter of fact, this is only true if one disregards the “trivial” part = {(p, q) ; p = k or q = k), which is 
the union of two ^f-dimensional hyperplanes. This “trivial” part gives a vanishing direct contribution to the collision 
integral (25i because either n(p) = n(k),n(q) = n{() or n(p) = n{(),n{q) - n(k), and it is thus generally ignored. 
However, we shall see below that it may be of indirect importance because any intersection of the “non-trivial” part 
with leads to sets of critical points, generically of dimension d - and possible divergences. 


The condition for a critical point on the 4-wave resonant manifold is 

Vw(p) = Vw(q) = Vcu(f) 


(18) 


with the group velocity of all three wavevectors p, q, ( the same. Degeneracy depends upon the rank of the D x D 
Hessian matrix 


V » = I ^ ® ® Vto(e) -V ® Vw(f’) 

^ _V®Vw(/’) V ® Vw(q) - V ® Vw(f’) 


(19) 


the critical point being non-degenerate if this matrix has full rank and otherwise degenerate. 

It is worth discussing the case of a general isotropic dispersion law a)(k), as a preliminary to some specific examples 
below. The condition for a critical point 

aj'ip)p = co'{q)q = co'({)( 


can be met in one of two ways: 

(i) \<jj'{p)\ — \(x)'{q)\ - W{()\ + 0 and the vectors p,q, ? are all collinear with k (parallel or anti-parallel 

depending on the sign of w'), 


(ii) (Jj'ip) - co'iq) = = 0 and no restriction on the vectors p, q, {. 


We shall refer to the first case as a non-null critical point with non-vanishing group velocity of the waves and to the 
second as a null critical point with zero group velocities. We include for completeness the third case 

(Hi) At least one of oj'ip), oj'{q), <jj'{i) is infinite. 


This is what we earlier termed a pseudo-critical point. These present generally no difficulty since the density of the 
phase measure with respect to surface area on the resonant manifold vanishes at pseudo-critical points, with an infinite 
gradient. Note also for the isotropic dispersion that 

V ® Vw(k) = w"(k)k ® k -H - k ® k), 

k 

from which it is easy to determine the rank of the Hessian matrix ( [T9] l. 

We now consider several concrete examples: 


2.2.1. Surface gravity-capillary waves 

An illustrative example is surface gravity-capillary waves with dispersion relation co{k) - sjgk ■¥ crk^, where g is 
the acceleration due to gravity and cr - S Ip, where S is surface tension. It turns out that most of the relevant features 
appear already for the idealized one-dimensional case d — 

Consider first pure surface gravity waves with dispersion law oj - This case has no null critical points. 

Because any non-null critical points must have wavevectors all parallel to k, we may check for their existence in the 
simplest case = 1. As a matter of fact, the resonant manifold is analytically known and explicitly parameterized 
for d - and consists of points for which one of the wavenumbers out of p,q,k,€ has an opposite sign from the 
others d. There are thus also no non-null critical points ind - \ and therefore none for > 1. There are, however, 
pseudo-critical points in = 1 where the nontrivial portion of the resonant manifold intersects the trivial parts, as seen 
in Fig. 7 below. These occur at the points where either p = 0 or ^ = 0 and u' diverges. The non-trivial part of the 
resonant manifold is a union of three smooth pieces, joined at the pseudo-critical points, but the phase measure on it 
is locally finite (and the 4-wave interaction coefficient zero Il22l '). 
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Figure 7: Resonant manifold ‘R^^\ k = \ for surface gravity waves in d = 1, with the non-trivial part plotted in gray and the trivial part in dashed 
black. 


However, if the surface tension effect is included, there will be critical points as well as pseudo-critical points, as 
seen in Fig. 8, due to the inflection point of u(k) - K null critical points are not present. 

However, for k k^ there is a distinct wavenumber { k satisfying a)'(k) - a>'{£), implying that there are two non-null 
critical points for p = k, q = f and q = k, p = f with t - fk. To determine their degeneracy, we again examine the 
Hessian ( [T^ . At least one of the diagonal term vanishes and the off-diagonal term -V ® Vw(/’) has eigenvalue -u)"{€) 
with multiplicity 1 and with multiplicity d - 1. If k K such that u)"{k) + 0, then + 0 and there 

exist two non-degenerate critical points; \fk - k, such that cj"ik) - 0, then cj"{€) - 0 and there exists one degenerate 
critical point. In Fig. 8, we show for d - I the typical resonant manifold for k < K (left panel), k - K (middle 
panel), and k > K (right panel). In general, the degeneracy 5 = 0 for k k, and 5=1 for k - k^. As we shall discuss 
below, this implies that the phase measure remains hnite for the physically relevant case d -2. 



Figure 8: The figures show the resonant manifold for surface gravity waves in d = 1 for g = cr, with the non-trivial part plotted in gray and the 
trivial part in dashed black. Here k = 0.3 for the left panel; k = J(2- V3)/ V3 for the middle panel; k = 0.5 for the right panel. 


It should be noted that the wavenumber k, ^ 0.4(g/cr)F2 expected to lie in the transition range between two 
energy cascades, one at low wavenumbers driven by quartet resonances of gravity waves and another at high wavenum¬ 
bers driven by triplet resonances of capillary waves muais]. The quartet at the critical points in Fig. 8 has two 
wavevectors with each magnitude k,{, and one of the magnitudes k,{ is always < k^ and the other > k^. Thus, any 
possible observable effect of the critical points would presumably be found in the transition region where k k,, 
perhaps in laboratory experiments where there is no large scale-separation between the gravity and capillary wave 
regimes. However, hniteness of the phase measure for c/ = 2 makes it unlikely that there are any appreciable effects. 
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2.2.2. Wave propagation along an optical fiber 

A very similar situation to the previous one, but for = 1, occurs for optical wave propagation along a fiber, 
modeled by a Id nonlinear Schrodinger (NLS) equation with third-order dispersion Il24ll25]l . The dispersion law i^ 

toik) ^ sk^ + ak\ (20) 


for s = +1, with an inflection point at k^ - -sl3a. The collision integral is 

Ck[n] - - ^dp dqd€ 5{p + q - k - {)6{ijj(p) + cj{q) - cj(k) - (x){()) 

Xn{k)n({)n(p)n(q) 


1 


1 


1 1 
n(k) n{€) n(p) n(q) 


( 21 ) 


It can be easily shown here that 

E*'‘^Hp, q\ k) = a)(p) + aj(q) - a>{k) - u{p + q - k) - 3a(p + q - 2kt)(p - k)(q - k). 

The non-trivial part of the resonant manifold is the straight line p + q - 2kt in the pq'-plane, independent of k, 
and there are non-null critical points at the intersections with the trivial part, at {p, q) - (k, 2kt - k) and (2A:, - k, k). 
The collision integral on the non-trivial part becomes 


Ck[n] 


1 r n{k)n{2k, - k)n{p)n{2k, - p) 

3^j \p-k\ \p + k-2K\ 

1 1 1 1 

X - + - 

n{k) n{2kt - k) n{p) n(2kt - p) 


( 22 ) 


As expected for D — 2d - 2, these critical points produce logarithmic divergences in the phase measure at the points 
p - k,2kt - k. For the special choice k - k^, there is a double pole at p = k^ when V ® becomes identically 
zero at the critical point (p, q) - {k,, k,). This degenerate critical point corresponds to a triple intersection between all 
three smooth pieces of the resonant manifold; see Fig. 9 below. 



Figure 9: Resonant manifold k = \ of one-dimensional 3rd-order dispersive optical waves for s = -1, cr = 3, so that k, = 1, with the 

non-trivial part of the manifold plotted in gray and the trivial part in dashed black. The triple intersection is a degenerate critical point. 


Because the critical points correspond to intersections with the trivial part of the resonant manifold, the integrand 
of the collision integral vanishes at those points and thus the integrals may be finite. For example, assume that n{p) 


^As usual in application of the NLS equation to optics, the space variable z and the time variable t have their roles exchanged, and thus also the 
roles of wavenumber k and frequency to. However, here we revert to the notations used elsewhere in our paper. 
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is twice-differentiable in the vicinity of k and let 6phe the scale of variation of n(p) near that point. Then one easily 
finds by Taylor expansion that the contribution to Ck{n\ from integrating over \p - k\ < bp is 




1 

67r|a| 


/ n"{k)rP'{2k^ - k) 

A 2 


{n'{k)fn\2k, - k) 
n(k) 


+ {k <r^ 2kt - k) 


{Spf 

\k-kA 


(23) 


with a similar contribution coming from p ^ 2k^ - k. When \k - kj\ < dp this is replaced by a contribution from the 
double pole: 

^ [n''{K)n\K) - 2{n'(K)fn\K)\ (6p). (24) 

These are finite as long as n(p) is twice continuously differentiable near the poles. More generally, the contribution to 
the collision integral is hnite if n{p) has cusp-like singularities near p = k,, n{p) - n{kt) ~ A\p - k,f with c > 1. 

Although the collision integral remains finite under the assumptions stated above, it may nevertheless become 
large, in the sense that the nonlinear frequency Tk{n) - Ck{n)ln{k) could be of the same order as (or greater than) 
the linear frequency toik). If so, this violates the condition Ykiri) cx){k) required for validity of the kinetic equation 
ifTSlfTOl l^. Indeed, it was shown in the numerical study ll24l that, for certain values of the cubic coefficient a, the 
ratio Tk{n)lu{k) exceeded 1 near k - k, and, in that case, there was no longer quantitative agreement between the 
predictions of the kinetic equation and ensemble-averaged solutions of the NTS equation (see Figs.2 and 5 in 124)). 
This is a physically interesting example which shows that resonance Van Hove singularities can lead to a breakdown 
in validity of kinetic theory, even when the collision integral remains finite. 


2.2.3. Electrons and holes in graphene 

We have so far considered classical 4-wave systems, but resonance Van Hove singularities can also occur in 
quantum wave kinetics. The collision integral of the quantum kinetic equation has the typical form 

Ck.s[n] ~ ^ f ^‘‘P ^‘‘‘1 -H f - p - q)b(Wi(k) + Wj'(0 - w,5«(p) - Wj'"(q)) 

x<;q'"'[(l ± «.(k))(l ± nA{))nAp)nMq) 

-nAk)n,,(()(l ± n.5«(p))(l + n.5'"(q))] (25) 

with + for Bose/Fermi particles, respectively. (E.g. see section 2.1.6 of 121). The presence or not of resonance Van 
Hove singularities in the quantum case is thus governed by the same considerations as for classical wave kinetics. 

A concrete example of physical interest is the dispersion law a)s{k) - csk, i = +1 which for d = 2 describes the 
band energies of electron-hole excitations in graphene near the Dirac points. The quantum kinetic equation has been 
used to predict electron transport properties of pure, undoped samples of graphene mini ED. The resonant manifold 
is a three-dimensional surface embedded in the Euclidean space of dimension D - 2d - A. Any critical points are 
clearly non-null with all wavevectors collinear. Because a>''(k) = 0, the Hessian matrix ( [T9| ) has null eigenvectors 
(k, 0), (0, k) and is co-rank at least 2. Thus the critical points are all doubly degenerate. In fact, the critical points lie 
on a 2-dimensional surface. Eor example, if k = (k, 0) with k > 0, then the critical set for s, s', s", s'" = 1 consists of 
p = (pi,0), q = iqi,0) satisfying pi > 0,q\ > 0, and pi + qi > k. 

In Eig. 10, we plot the resonant manifold "RJA A k = ei for 4-wave interactions in graphene. Eor visualization, we 
employ the notation in Il27l[1411151 that the four resonant wave vectors k, t, p, q are instead denoted by k, p, k+q, p-q. 
The resonance condition in these variables becomes |q -H k| -H |q - p| - k + p, which, as noted in Il27l . corresponds at 
fixed p values to ellipses in the q-planes, with foci at -k and -i-p. When k||p, these ellipses degenerate to line-segments 
[-k, H-p], whose set union comprises the critical set. In Eig. 10 we plot 2D sections of the resonant manifold in the 
3D space (q^, qy, Px) at three different values of py. Eor the chosen value of k = Ci, the 2D critical set is located on the 
Py - 0 section, as shown in the middle panel. 

Eor graphene with D = 4, this 2-dimensional critical surface produces a logarithmic singularity in the phase 
measure, as has been previously noted ll27l[T4lfT5l . Eor example, for k = kei,k > 0, using the notation k, p, k-i-q, p-q 
for the wavenumber quartet, and writing p = ipi,p±), q = (qi,q±), it was pointed out in m that to quadratic order 
in the transverse variables near the critical set 

A _ ql _ (p± - q±f 

2p 2(k + p) 2(p - q) 
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kH-p- Iq-i-kl - |q-p| 












Figure 10: The section of the resonant manifold k = (1,0) for fixed py (py < 0 for the left panel; = 0 for the middle panel; p_y > 0 

for the right panel), with the non-trivial part plotted in gray and the trivial part in dashed black. The 2D critical set in the middle panel is plotted in 
dark gray. For py 9 ^ 0 the horizontal sections of the resonant manifold at fixed p.v-values are ellipses, but for py = 0 and px > 0 these sections are 
line-segments. 


" ~ 'll (26) 

2(k + q)(p - q) 

where (12 are the roots of the quadratic polynomial in q_i_lp_i_ dehnedby the first line. At each fixed value of pi,qi that 
corresponds to points in the critical set, the integral over the transverse variables p_L,qi_ is logarithmically divergent. 
This can be seen most easily by changing variables to 


rj_ ^ qj_-( iPj_, Sj_^q_L-(2P± 


with the Jacobian of transformation 




d(q±,p±) 


1^1 -^2l 


k(k + q)(p - q) 


k + p 


Thus, 




dq^. + p - |q + k| - |q - p|) 


p{k + q)(p - q) 


1*^1 


dsx_ d(rxij_) 


and the latter integral exhibits, for each hxed value of pi,q\ corresponding to points in the critical set, the same type 
of logarithmic divergence that was observed for Rossby/drift waves 'md -2. However, electron-hole interactions do 
not vanish near the critical set, unlike the case for Rossby/drift waves, so that the divergence is unregulated. 

The local logarithmic divergence can be removed by resonance-broadening of the delta function (e.g. |[3l, section 
6.5.2). Assuming for simplicity a constant resonance width y, the delta function (5(rxSj_) is replaced by the Lorentzian 
(Cauchy distribution) 

y 1 

-r. 

71 

Since dyCr^s^) ~ l/Try for « y, an integral over a small neighborhood of = 0 is now hnite. The 

quantitative behavior can be seen from the following integral 

C* 

/ = ' 


r“ r r dri. 

I dr 2. I ds2 dy(rxij_) = 2 j - 


jj_r2 r 
Jo J-b 


ds± 


b 4 + (r/rj2 


with cutoffs fl, b. The region of integration ylb gives a y-independent contribution, whereas for values 'S> yjb 
the inner integral is ^ 1 and thus 

I ~ 2 ln(ah/y), y ab. 

This asymptotic evalation can easily be made rigorous by noting that I - ^Ti 2 (ah/y) where Ti 2 (x) = du “ 
is the inverse tangent integral (HS), Ch.VII, §1.2). The cutoffs a,b may come from limits to the magnitude of all 
wavevectors, e.g. the size of the Brillouin zone in graphene. The cutoffs a, b may also be associated to the maxi¬ 
mum size of the region where ^(p, q; k) oc We shall discuss in further detail in section |^how the logarithmic 
singularity is understood to affect the electron-hole kinetics in graphene. 
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2.3. Summary 

This section has shown by various concrete examples that critical points in the resonance condition (resonance 
Van Hove singularities) occur in many common wave kinetic equations. They lead to geometric singularities in 
the “resonant manifold”, which is thus no longer a true manifold. The singularities may furthermore lead to local 
non-finiteness of the phase measure appearing in the collision integral, which is associated physically to the infinite 
scattering time for locally non-dispersive waves. Such a diverging phase measure may nevertheless produce a finite 
collision integral, e.g. due to a vanishing interaction coefficient (Rossby/drift waves, d - 2) or due to cancellations 
between terms in the collision integrand (waves in an optical fiber, d = 1). When the collision integral itself diverges 
or even if it is finite but large, standard kinetic theory may break down (optical waves, = 1; electron-holes in 
graphene, d - 2). Before we discuss this latter situation in section]^ we first discuss in the following section more 
generally the conditions under which a critical point leads to a locally infinite phase measure at the singularity. 


3. Phase Measures and Their Finiteness 


While the previous section considered a rather disparate set of examples, the present section develops some quite 
general results about the effect of critical points on the local finiteness of phase measures. For the case of non¬ 
degenerate critical points we can carry over essentially unchanged the considerations of Van Hove in his classic study 
El. Thereafter we discuss briefly the case of degenerate critical points. 

We begin with a mathematical issue that we have neglected until now: the resonance function ^(p; k) is generally 
not smooth enough in the phase-space variable p = (pi,..., p^_ 2 ) in order to make the delta-function i5(£’(p;k)) 
meaningful in any naive sense. Thus, the phase measure that we have defined formally by dfi — d^p (5(£(p; k)) has 
no actual mathematical meaning. To give it a proper definition, one must return to the derivation of the wave kinetic 
equation. A standard multiscale perturbation argument Ema shows that what appears in the kinetic equation is 
actually an approximate delta function of the form 

dH£(p; k)) - ^sinc2 (|£(p; k)) , (27) 

for a time T chosen so that w(k) » l/T » rk(n), where “sine” denotes the cardinal sine function sinc(.r) = smxjx. 
A physically motivated definition of the phase measure is thus as a suitable limit 

dp - lim d’^p k)). (28) 

T —*oo — — 


We can show using a Daniell integral method that this limit yields a well-defined measure which is unique among 
those absolutely continuous with respect to D-dimensional Hausdorff measure TT® ' on the resonant manifold ‘Rk and 
that this measure satisfies 


t/'K^-'(p) 

|V£(p;k)r 


(29) 


Because the details are rather technical, we provide them in another paper If30l . An alternate approach is suggested 
by field-theoretic derivations of wave kinetics Gl, which instead yield an approximate delta function of Lorentzian 
(Cauchy) form: 


(5y(£'(p;k)) 


y 1 

7T ^^(p; k) H- ’ 


(30) 


and which suggests to take a similar limit y ^ 0. This was the starting point of Lukkarinen & Spohn ll3Tll to define 
the phase measure, in a slightly different context. We also consider this alternate approach in Il30l and compare with 
the method of OTlI . 

Now using the relation ( [29] l we study the local finiteness of the phase measure p in the neighborhood of a non¬ 
degenerate critical point, following the basic idea of Van Hove El. who analyzed this question for the energy density 
of states using the Morse Lemma. Since k appears simply as a parameter in the argument, we omit it and write simply 
^(p). We recall here the statement of the Morse Lemma IIIIIOI: if ^(p) has a non-degenerate critical point p such 
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that E is k>\ with respect to p in a neighborhood of p , then there is a diffeomorphism ^ of a neighborhood 
t/ of p with a neighborhood V of 0 = ^(p ) such that £(q) = ^(^^^(q)) has the canonical form 

£(q) = £(pp - Iq’l" + Iql' 

where (q”,q^) = q = ¥>(?) with q"^ e K®* and D + D+ - D. It is possible that D - 0, in which case £(q) = 
£(P ) + IqP or that D+ - 0, in which case E{q) - E(p ) - |q|^. Note that D± are just the number of positive/negative 
eigenvalues of the Hessian of £ at p . Van Hove assumed further that ip can be chosen to preserve D-dimensonal 
volume (Lebesgue measure), in which case 

/= r <5(£(p)) = r d’^q5{E{q)), 

Ju ~ Jv ~ 

and the righthand side can then be shown to be finite/infinite by a direct calculation. It has indeed been proved 
subsequently that such a volume-preserving (“isochoric”) choice of (p is possible, if £ is a C“ function l32l . However, 
even if £ is only the neighborhoods U, V can always be chosen so that the Jacobian determinant of p satisfies 

0 < c < |det(D^)(p))| <C < oo, pet/ 

and since 

r d’^qdiEiq))^ f d^pldet(Dp)(p))l6(E(p)) 

Jv ~ ~ Ju ~ ~ ~ 

the two integrals J = d(E(q)), I - j^d^p 5(£(p)) are either both finite or both infinite. Thus, the question 

again reduces to an elementary calculation of the integral J. 

Taking therefore a critical point p on the resonant manifold, satisfying £(p ) = 0, the condition for the resonant 
manifold in the q-coordinates in V becomes simply 


Iq'l 


= Iq1. 


When either = 0 or D+ = 0, the resonant manifold reduces to the isolated point q = 0, and one can easily show 
that 7 = 0. Thus we assume D± + 0. We use 


7 = 


r A5(£(q))= f 

Jv ~ JvnK 


7’K"-i(q) 

|V£(q)| 


and V£(q) = 2(-q ,q+), so that |V£(q)| = 2 V2|q | on !R. To simplify the calculation, without loss of generality, 
we take the neighborhood V to be a Cartesian product of two balls of radius q, V - B{(r,q) x B(Q^,q). Using 
hyperspherical coordinates for both q“ and q^, the calculation reduces to|^ 

7 = r d'H^ ^tf'} — 1=— r 777^+"'(q^) = (const.) f q’^^^dq-, 

JB{0 , 7 ]) ~ 2 \2.Cl- Jq+~q- ~ Jo 

where a factor oc * arises from the (D+ - l)-dimensional Hausdorff measure of the q^-sphere of radius q_. Clearly, 
7 < oo if D > 2 and 7 = oo if D < 2, exactly as had been concluded by Van Hove for the energy density of states. 


^ The Fubini-like theorem that we use here follows from the co-area formula of geometric measure theory. E.g. see Theorem 3.2.22, ES 
For any W c which is M-rectifiable and ‘H^-measurable, Z c R^, N < M which is A/^-rectifiable and Ti^-measurable, a Lipschitz map 
f : W ^ Z, and any non-negative, -measurable function g : VF —» R, J^g{w)(Jf)(w)d']-(^(w) = g{^')dd-(^~^(w'). 

We apply that theorem with M = D - I, N = />_, W = V 0 V = Z x Y, Z = B{0~,r}), Y = / : W —> Z is the restriction 

to W of ;r : y —^ Z, the projection onto the first factor of V, and g = 1/|V.E|. Note that for each q“ 6 Z, the q"^- 

sphere of radius q- centered at 0“^, and that the Jacobian Jf = 1. Finally, W is obviously (D - l)-rectifiable and Z is Z)_-rectifiable. Thus, 
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In particular, we find a logarithmic divergence for D - 2, when one can write simply q = (q^, q+). Introducing the 
new coordinates 


r^q^-q_, s ^ q+- q- 


with Jacobian of transformation I , = 2 one can write 

\d(q+,q-)\ 


E - q\ - q^_ - rs 


and thus 


7 = 


/ 


dq+dq- 6{q_^ 


i) = i 


f 


dr ds 6(rs). 


I\q\<7j ^ Jr^+s^Klrf 

We see from this formula that our previous considerations on the effect of resonance broadening for graphene carry 
over to the general case, with the local divergence removed and replaced by a logarithmically large value cx ln(af>/'y) 
for resonance width y. 

We have thus obtained quite general results on the local finiteness of the phase measure in the vicinity of a non¬ 
degenerate critical point. These general considerations explain the specific results we found in earlier examples, such 
as the logarithmically divergent phase measure for 3-wave resonance of Rossby/drift waves in 7 = 2 (D = 2) and the 
finite phase measure for 4-wave resonance of capillary-gravity waves ind - 2{D - 4). 

Let us now consider briefly the effect of degeneracy. The simplest situation is when the critical points lie on a 6- 
dimensional submanifold where V£’(p) = 0, which immediately implies a degeneracy degree of at least 6 at each such 
critical point. This corresponds to the situation where in a neighborhood U of each critical point p , with ^(p ) = 0, 
there is a diffeomorphism ip with a neighborhood V of 0 = ^(p ) such that 


£(q) 


P + Iq1" 


where (q°, q , q^) = q = ^(p) with q* e M®* and q° e K"*. In this case, D+ + - D - 6 = D'. One can take the 


neighborhood V, without loss of generality, to be of the form V - Vq x V', where Vo is a neighborhood of 0° e 
V' is a neighborhood of (0^, 0^) 6 In this case, it is seen that 


and 


J = £^(Vo)J', f = 


'= r d’^'q 

Jv ~ 


5(£(q')) 


with q' = (q“, q^) and the d-dimensional Lebesgue measure (volume). Here J' has the same form as did J for the 
non-degenerate case, but with D’ replacing D. Thus, the measure is locally finite at each critical point for D' > 2 but 
locally infinite for D' < 2. 

This analysis explains the results we obtained in several concrete examples, such as the logarithmically divergent 
phase measures for 3-wave resonance of inertial waves in 7 = 3 (Z) = 3, d = 1) and for 4-wave resonance of electron- 
hole excitations of graphene m d - 2 {D - A, 6 - 2), both with D' - 2. 

Another case of interest is an isolated critical point with degeneracy degree 6 (co-rank of the Hessian matrix). 
The classification of isolated critical points for differentiable functions E belongs to the field of singularity theory; see 
El- We shall not discuss such a classification in detail here, but we briefly mention in this light the doubly degenerate 
critical point obtained for wave propagation along an optical fiber, from section 2.2.2 For the distinguished value 
k — K (known as the “zero-dispersion frequency” in the nonlinear optics community) one can express E{p,q-,K) 
locally at the degenerate critical point (k*, k*) in terms of the deviation variables 5p - p - K, dq - q - K as 

E{6p, 6q', K) - 3a(6p + 6q) 6p 6q. 

After a linear transformation 6p - (x- y)f(6ay^^, 6q - —{x -i- y)/(6a)''^^ this becomes 

E(x,y;k^) - x^y-y^. 

The latter is the normal form for the family D7 in the classification of “simple” singularities for differentiable real 


functions in 113411 . The general possibilities are quite rich and complex, and still the subject of mathematical inves¬ 
tigation. However, we note from the optics example that the effect of degeneracy is again to worsen the divergence 
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of the phase measure, whose density now exhibits a double pole rather than the simple poles (leading to logarithmic 
divergences) found for the non-degenerate critical points when k + k^. 

The cases of degenerate points that we have discussed here are by no means exhaustive. For example, one could 
have a set of critical points with degeneracy degree 6 comprising a Do-dimensional submanifold with 1 < Do < 6. 
There is also the possibility of pseudo-critical points, but, as discussed earlier, these will produce no divergence of 
phase measure unless the geometric singularity is so severe that the resonant manifold develops a locally infinite 
Hausdorff measure. 

The general morals to be drawn from our discussion are as follows. Pseudo-critical points should usually yield 
a locally hnite phase measure and are “harmless” for kinetic theory. True critical points are potentially “dangerous” 
and can lead to locally infinite phase measure, especially in situations of low dimensions d, low-order N of resonance, 
and/or high degeneracy degree 5. The divergence of phase measure due to such singularities can be rendered harmless 
by vanishing interaction coefficients or by cancellations in the collision integral. In the next section we explore the 
opposite situation when the standard collision integral with exact resonances remains divergent. 

4. Singular Wave Kinetics 

We have seen several examples (inertial waves in c/ = 3 fluids, optical wave propagation along a = 1 fiber, and 
Dirac electron-hole excitations in d - 2 graphene) where an “unprotected” resonance Van Hove singularity leads to 
a breakdown of standard wave kinetics. This is analogous to the situation in the theory of low-amplitude acoustic 
wave turbulence, except that now the breakdown of dispersivity of the waves occurs only locally on the resonant 
manifold rather than for all resonances. In the case of acoustic waves, generalized equations were derived to describe 
the “singular wave kinetics” of semi-dispersive waves, both by multiple time-scale perturbation theory IfT^ and by 
field-theoretic methods lua. One can expect that such singular kinetic theories will apply more generally, even when 
the critical set is only a subset of the entire resonant manifold. We shall briefly illustrate this situation with the example 
of electron-hole kinetics in graphene. 

We first discuss graphene from the point of view of multiple time-scale perturbation theory, as originally devel¬ 
oped in lfT2l . It should be stressed at the outset that the applicability of perturbation theory is itself a nontrivial result, 
because electron-hole excitations in graphene on most substrates are not weakly coupled (nor inflnitely-strongly cou¬ 
pled). It is instead believed that the coupling becomes weak for sufficiently low wavenumbers because of a globally 
attractive, asymptotically-free renormalization-group (RG) fixed point. See llTSl [3^iJTll for discussions of this issue. 
Thus, the perturbation theory argument must be applied to a low-wavenumber renormalized theory. A further very 
significant complication is that the Coloumb interaction is long-range and many-electron effects such as dynamical 
screening are expected in graphene Il38ll . and these effects do not appear at any finite order in naive weak-coupling 
perturbation theory. Although such a naive perturbation theory analysis is therefore very incomplete, we And that it 
provides useful insight. 

The perturbative derivation of the quantum Boltzmann equation is much the same as the derivation of the kinetic 
equation for weakly-coupled classical waves, e.g. see ||6]|2l. Start with a general quantum Hamiltonian H - Hq + Hi 
with free part 

= Z f ‘^kajAk)rl(k)j,,Jk) (31) 

s,a ^ 

and interaction part 

" I X X r ^‘^kl^^2C^‘^k3d^k4 5‘^(ki +k2-k3-k4) 

5i,5'2,-s'3,-y4 ^ 

X 7'k;£rk.<.(k4)rL(k3)r..a(k2)r../.(ki), (32) 

where ym(k), 7 jQ(k) are standard creation/annihilation operators. We assume here that these operators obey canonical 
anti-commutation relations, as appropriate for fermionic electron & hole excitations. The indices a, b are summed 
over integers 1,..., A, to incorporate a possible A-fold degeneracy of the fermions. For an infinite-volume system. 
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define the mean occupation number ni(k, f) by (710(14, ^)7.!'a'(k^ 0) = ^^(k, t)5ss’6aa’ 5'^(k - k')- Perturbation theory in 
the small parameter a yields the result that 


n*. (ki,f)-nj,(ki,0) 


■^4 


J d^k2d^k3d^k4 5''(ki+k2 


•k4) 


:^51^2'S3-^4 ^2 


''^kik2k3k4'' 

k,[”"i (*^1 )”'2(k2)ni3(k3)nj,(k4) - n^, (kj (k2)n,3(k3)nj,(kr)] 


where 


pS\ i'2^3-54 

'^kik2k3k4 


- uji,(ki) + Wi2(k2) - (Us3(k3) - Wj^(k4), 


1 , 


pSiS2S2S4 _ ^.IT 

^kik2k3k4 “ 2 ' 


51 52i'3-S4 

ki k2k3k4 


-T, 


51 52 54-^3 [2 
kik2k4k3 1 


A^- 1, 


rS\S2S2S4 [2 

'kik2k3k4l 


A^- 1, 


pj'l52-54‘S3 [2 
-kik2k4k3' • 


and 


n3(k, f) = 1 - n^Ck, t). 


(33) 

(34) 

(35) 

(36) 


See Appendix A for the details of the derivation. An elementary calculation gives 

1 - cos(Et) 


£2 


= nt 6,(E) 


where dt(E) is the approximate delta function in eq.(|27|i. When the condition ^kikjkski “ ® defines a non-degenerate 


resonance, the righthand side of (33 i exhibits a secular growth cx a^t. This secular behavior is removed by choosing 
the occupation number to satisfy the quantum kinetic equation on the slow time scale t = a^t, with collision integral 

This standard derivation fails for electron-hole kinetics in graphen^ when d -2 and 

ii>s(k) — svfk, s = +1. 

The index A = 4, for the two electron spins and two Dirac points (valleys) in the Brillouin zone. For the expression 
^k|k2k3k4 rising from Coloumb interaction, see eq.(3.9) of ifTsl . The degeneracy of the condition ® 

changes the asympotics of the integral in eq.(33i. Note that momentum and energy conservation allow non-trivial 


resonances only for electron-electron/hole-hole collisions (all i’s of the same sign) or electron-hole collisions (one 
i = -1-1 and one s - -1 in both incoming and outgoing states). For simplicity, we discuss explicitly here only the 
first case. To obtain the long-time asymptotics, it is useful to take two derivatives with respect to time, to obtain the 
contribution with all i, = s: 


^[nAk,t)-nAk,0)]^-2a^ 


I 


d% d% cos(£4Y^;^f)M++++(k3, k4), 


where k 2 = k 3 -i- k 4 - k and 


(37) 


M"+"+(k3,k4)_=£^,2kA4 

X [nj(k)n,(k2)n3(k3)n3(k4) - nj(k)n,(k2)n3(k3)nj(k4)], (38) 


P+ +++ 

^kk2k3k4 


_ |y^ ++++ 

“ 2 ^ ^^2k3k4 


+ |2 

^ kk 2 k 4 k 31 


+ (A-l)|7’, 


kk2k3k4 


(39) 


^The physical dimensions require a bit of discussion. The 1-particle energies in the free part of the Hamiltonian are ^’^(k) = /i<£» 5 (k), so our 
units correspond to h = 1. The Coulomb potential is y(r) = jKr, where k is the background dielectric constant. Hence, the parameter a in the 
interaction Hamiltonian for graphene corresponds to the dimensionless “fine-structure constant” a = jKhvf, in units where v/r = 1. We have kept 
the factor v/r in the dispersion relation, although it is strictly just 1 in our units. 
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Using polar coordinates k = k(cos 6, sin 6) with fixed direction angle 0 of wavevector k and k, = ki(cos(6 + Oi), sin(6> + 
9i)), i = 3,4, this can be rewritten as 


■^[risiK t)-ns{K 0)] 


/^OO ^oo ^JT r^TX 

I ki dkj, I ^4 dk4 I dO^ I d04 
*Jo «Jo *J -n k)-TT 


g'f'^k,k3+k4-k.k3.k4^++++(;t3, k4, 03, 04), 


where a simple calculation gives 


k,k3+k4-k,k3,k4 - ^3 + ^4 k 


- slikj + k/[- k)^ - Ik^k^il - cos(03 + 04)) + 2^3^(1 - cos 03) + 2 k 4 k(l - cos 04) 

0^H0 


to quadratic order in 0 = (03,04), with the 2x2 matrix 


H = 


1 


^3(^4 — k) ^3^4 


2{k^ + k4- k)\ ^3^4 k4(ki — k)j 


kk-ikd 


(40) 


(41) 


Using det//- 4 (k-k^ 
of stationary phase 1391 . giving 


, the integral over angles in eq.(40l is then evaluated asymptotically for t ■ 


(42) 

I by the method 


^[n4(k,f)-ni(k, 0)] 


In 

t 


If 


^3^4 dkj,dk 4 |det//| 02 M^^^^(03 = 04 = 0) 


In 
+ — 
t 


k-i >0,k4>0,k2 +^4 >k 


'fi 


^3^4 dk^dk4 I detH\ 02 M^^^^(03 = 04 = 0) 


_ _ 2 a^ 1 rr 

^ t JJk, 


:3>0,k4>0,ks+k4<k 

dkj, dk4 


+ 0{n) 


k^ >0,^4 >0,1^3 +1^4 >k 


(k^ + ^4 - k)kj,k4 


m++++( 03 = 04 = 0) + 6)(r0. 


This integrates to 


n,(k, t) - n^(k, 0) 

.2 




dk-i dk4 


3 >0,A:4 >0,^3+A:4 >A: 


{k^ + ^4 - k)kzk4 


M++++(03 = 04 = 0) 


(43) 


(44) 


as f ^ 00. The crucial point is that the leading secular growth is now faster than f by a logarithmic factor Inf, as a 
consequence of the degeneracy of resonance. Note that the entire contribution arises from the critical subset of the 
resonant manifold, with all quartet wavevectors parallel to k. 

Taking into account the electron-hole scattering changes this result only by appearance of an additional term, 
which is obtained by a very similar calculation. The complete asymptotics is given by 


n,(k, t) - n^k, 0) 
2op- 


^finf rr 

2!- JJk, 


dk-i dk4 


k-i>0,k4>0,k^+k4>k 


(/r3 kiX — k^k'},kiX 


M(03 = 04 = 0) 


(45) 


with 

M(k3, k4) = M++++(k3, k4) -h M+"""(k3, k4) 
and the electron-hole scattering contribution 

M-"ik3,k4) = <;t;k3,_k3 

X [n,(k)n_j(-k4)n,(k3)n_,,(-k2) - ns(k)n_j(-k4)nj(k3)n_,(-k2)], (46) 
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P+-+- 

^kk2k3k4 


I^42k3k4 


J+ —+ |2 

^ kk2k4k3 1 


+ (A^-1)|7’, 


kk2k3k4 


+ {N-\)\T, 


kk?k4k3 


(47) 


We have made a change of variables k2 -k4, k4 —» -k2 in the integral for the electron-hole contribution so 
that the range of integration is the same as for the electron-electron contribution. The physics of the electron-hole 
scattering term is easy to understand, if one recalls that a hole excitation with wavenumber k has a group velocity 
Va;_(k) = -v/rk which is opposite to the group velocity V<u+(k) = -HVfk for an electron excitation with the same 
wavenumber k. Hence, non-dispersive interactions with identical group velocities for a quartet of modes requires that 
the holes have wavenumbers anti-parallel to the wavenumbers for the electrons. The fact that electron-hole scattering 
couples occupation numbers for anti-parallel wavenumbers will be seen below to have interesting consequences. 

The leading secular growth in ( [45] ) can be removed, following the ideas in ifT^ . by allowing the occupation 
numbers to evolve on the time-scale ti = Inf according to the singular kinetic equation: 


dnAk)ldTi^-- ff dhdh / fe+^4 k)k^ 

>^{^kk^2k^|k4 [ni(k)ni(k2)n.v(k3)nj(k4) - ni(k)ns(k2)n.v(k3)nj(k4)] 

+^k,”k7,k2,-k2 [ns(k)n-.s(-k4)n,(k3)n_,(-k2) - ni(k)n_,(-k4)nj(k3)n,,(-k2)] |, (48) 


where in the collision integral all wavenumbers k, = A:,k are parallel to k and k 2 - k^, + kti - k. Similarly as in the 
work of Newell & Aucoin ifl^ on acoustic turbulence, this new kinetic equation is actually a continuum of uncoupled 
equations, one for each line along +k. Whereas the usual quantum kinetic equation holds on a time-scale T2 = a^t, 
the singular kinetic equation is valid at logarithmically shorter times. Put another way, for ti = Oil), T 2 = (9(1) 
one finds f Ti/a^ln(l/a) « T2/a^, when ln(l/a) » 1. It is important to note that the singular kinetic equation 
( |48| ) coincides (up to a constant of proportionality) with the result previously derived for electron-hole kinetics in 
graphene by mna, who made a leading-logarithm approximation to the divergent collision integral in the standard 
quantum Boltzmann equation. However, such a “derivation” of ( [48) l is inconsistent, taken literally, because it employs 
the quantum Boltzmann equation in a regime outside its validity. We discuss further below the derivation of nma, 
which must be consistently understood within a proper field-theoretic framework. Both our derivation and that of 
IHEl have also made an ad hoc assumption that the Coulomb interaction is dynamically screened at very low 
waveumbers, in order to eliminate an infrared divergence of the collision integral R but a proper derivation of this 
effect requires a more sophisticated many-body theory. 

At a time f = ti /a^ ln( 1 /a) with r i » 1 (but with t <Si T2 /a^), the solutions of the singular kinetic equation should 
be expected to approach a local equilibrium separately along each line in directions +k. This is what occurs in the 
singular kinetics for acoustic turbulence IfT^ and was also argued to occur for electron-hole kinetics in graphene in 
nma. The local equilibria of (|48|) can be easily checked to be of a generalized Fermi-Dirac form 


n,(k) = 


1 


exp[(y8j(k)ik H- fii(k))] -H 1 


s - ±1, 


(49) 


where 

^,(k) = /3e(k) -H j0o(k)i, /i,(k) = fi(k) + 4f(k)i 

are distinct inverse temperatures and chemical potentials, independently specified for each direction k and for electrons 
(s = 1) and holes (s = -1), subject to the conditions that jSe(k) must be even and ySo(k) odd: 


;6,(-k)=/3,(k), /3o(-k) ^-/3oik). 


(50) 


This last restriction arises from the vanishing of the electron-hole contribution to the collision integral in (481, whereas 
the same local Fermi-Dirac distribution ( |49| ) causes the electron-electron/hole-hole term to vanish without any restric¬ 
tion on parameters. Note that y6o(k) is the temperature asymmetry between electrons and holes, and 4f(k) is the 


^This divergence due to unscreened Coulomb interaction is simply exhibited in the coordinates used in Fig. 10, with ki = k, k 2 = p, k 3 = k + q, 
kt = p - q, for which the singular collision integral is with respect to the measure dpclq over the range p > 0 and p > q> —k. Since 

the d = 2 Fourier transform of the d = 2 Coulomb potential is oc Ijq, the interaction coefficients oc l/q^, leading to an integral over 

dqjq^ divergent a\.q = Q 
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chemical potential asymmetry. If one confines attention to solutions satisfying particle-hole symmetry appropriate to 
zero doping, n-s(k) = nj(k), then the possible equilibria are reduced to 

1 

nj(k) = -7-7-, i = ±1. (51) 

exp[(;S,(k)^+;ir(k))i] + 1 

Nonzero values of symmetric chemical potential /^(k) or of temperature asymmetry /lo(k) explicitly break particle-hole 
symmetry. The results (|49|),(|50ll,(|5T]) do not seem to have been given earlier in the literature, although they are implicit 
in the papers m\iB 1261 on dissipative transport by electrons in graphene. These local equilibrium solutions are not 
only stationary solutions of ( |48] ) but in fact should be global attractors. It is straightforward to show by standard 
arguments that the entropy 


5(6, = -.. 2 X“ 


kdk [«j(k) In nj(k) -H nj(k) ln«j(k)] 


satisfies an //-theorem of the form dS{V)ldT\ > 0 for solutions of the singular kinetic equation (481, separately 
for each direction k, and the only solutions with vanishing entropy production, dS (k)/(/Ti = 0, are the generalized 
Fermi-Dirac distributions (|49ll,(|50ll. 

Until now our discussion has been quite parallel to the case of acoustic turbulence lfT2l . but we now encounter 


a difference, because the critical set which dominates the singular kinetic equation (481 is not the entire resonant 
manifold for electrons in graphene, unlike the situation for acoustic waves. Thus, there are additional secular terms oc 
a^t in eq.( 331 which arise from integration over the remainder of the resonant manifold. To remove those secularities, 
one should impose an additional dependence upon the time variable T2, equivalent to the condition that satisfies 
the ordinary quantum Boltzmann equation on time-scales T2 = 0{\), or that c/nj(k)/c/T2 = Ckj(n) with the collision 
integral in ( [25| . This collision integral diverges for general distributions «s(k), but we can show that it is finite for 
local Fermi-Dirac distributions of the form (j4^. (Details will be presented elsewhere IHbI .) A simple picture thus 
arises for electron kinetics in graphene as a three time-scale problem. At the shortest times f = tq of order the linear 
wave period, the dynamics is dominated by the free part of the Hamiltonian, whose dispersive character drives the 
system into a locally quasi-free state completely characterized by the occupation numbers nj(k). Cf. the discussion in 

& section 9. At times t - TxjcP- ln(l/a) » tq the occupation numbers evolve further into the local Fermi-Dirac form 
, completely characterized by the local thermodynamic parameters y6o(k),/?<,(£),//(k),;(f(k) along each wavevector 
direction. Finally, at times t - » Ti/cr^ ln(l/Q') the distribution further relaxes according to the standard 

quantum Boltzmann equation. Because the ratio of times t with T2 = (9(1) and t\ = (9(1) is only logarithmically large, 
it is possible that the occupation numbers will not fully relax to a local Fermi-Dirac form ( [49| for times T2 = (9(1) 
and there may be corrections of order l/ln(l/a'). If there is no external driving or boundary conditions to keep the 
system in a dissipative non-equilibrium state, the subsequent evolution by the standard quantum Boltzmann equation 
will relax the system at times T2 » 1 to a global Fermi-Dirac equilibrium 


n5(k) 


expifisk H- yu) -I- 1 ’ 


with uniform values of inverse temperature [5 and chemical potential p. 

So far we have discussed electron kinetics in graphene from the multi-time perturbation theory viewpoint devel¬ 
oped by Newell-Aucoin IfT^ to describe semi-dispersive acoustic turbulence. There is however another point of view 
on singular kinetics for acoustic waves which was developed by L’vov et al. 03 , based on a Martin-Siggia-Rose field- 
theoretic formulation. In this approach, the starting point is a set of Schwinger-Dyson integrodifferential equations, 
which are exact and non-perturbative but non-closed. By a set of rational approximations based on weak non-linearity 
and self-consistency, the authors of na showed that the Schwinger-Dyson equations for acoustic wave turbulence 
can be simplified to a generalized kinetic equation. This has a form similar to the standard 3-wave kinetic equation 


(9,«(k, f) 


367r 2 r d‘^k2 f d‘‘ki |//§2dr(k)(s ■ w(k)) S‘^(s ■ k) 

5 =(- 1 , 52 , 53 )"^ 

x{n(k2)n(k3) - S2n(k)n(k3) - S3n(k)«(k2)}. 
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but with (5r(k)(w) a resonance-broadened delta-function or Lorentzian of the form (301, where 


r(k) = r(k) + r(k2) + r(k3) 

is a triad-interaction decay rate. The individual rate -/(k) is given by the imaginary part of the self-energy function 
Z(k), so that it must be determined self-consistently in terms of the solution n(k) of the generalized kinetic equation 
and is (9(e^) in the nonlinear interaction strength e. See eq.(B13) of IT?! . The collision integral of this generalized 
kinetic equation remains finite and free of any divergences due to the Van Hove-type singularities in the “manifold” 
of exact resonances for acoustic waves. 

This field-theoretic point of view is very closely related to the previous derivations of the quantum Boltzmann 
equation for electron kinetics in graphene iiaiia, which were based on a nonequilibrium Schwinger-Keldysh held- 
theory approach (the quantum analogue of the Martin-Siggia-Rose held-theory for classical dynamics). In fact, the 
authors of nana assumed that self-energy corrections will cut off the divergence of the standard collision integral 
due to the resonance Van Hove singularity, but without assuming an explicit form for these corrections. It is likely 
that the self-consistent approach of ns can be carried over to electron kinetics in graphene, e.g. with a wavenumber- 
dependent broadening corresponding to a quartet-interaction time 


r(k) = y(k) -H 7 (k 2 ) -H ^(ks) H- y(k 4 ). 


As discussed in section 2.2.3 such a broadening should cure the logarithmic divergence, if the resonance width is 


non-zero in the vicinity of the critical set. There is considerable interest in investigating explicit self-energy regular¬ 
izations, since it has been estimated that corrections to the leading-logarithm approximation may make a 30% change 
to the electrical conductivity of defectless graphene at experimentally realizable temperatures lfT4l . In addition to the 
interest for potential electronic applications of graphene, such a study would also help to assess the validity of the¬ 
oretical approximations for acoustic turbulence which, to our knowledge, have never been subjected to experimental 
test. More generally, electron kinetics in graphene is a problem which should illuminate the subject of singular wave 
kinetics, with applications to a wide variety of systems. We are currently pursuing such investigations Ho). 


5. Conclusion 

While we have developed no comprehensive theory for existence of resonance Van Hove singularities, the exam¬ 
ples presented in this work indicate that they occur rather commonly. Their presence may be due to disparate causes, 
including periodicity of Fourier space, anisotropy of the wave dispersion relation, or intersection of the trivial and 
non-trivial parts of the resonant manifold for 4-wave resonances. The basic requirement for such critical points is 
that there be distinct wavevectors with the same group velocity, which is facilitated by dispersion laws with segments 
strictly linear in wavenumber or with inflection points. Our several examples have presumably not exhausted the 
possible mechanisms to produce such resonance singularities. 

The effects of the singularities on the kinetic theory can range from none at all, to moderate, to quite destructive. 
As a general rule of thumb, the singularities are more threatening in low dimensions (d = 1,2). For example, the non¬ 
degenerate critical points for the gravity-capillary wave system in = 1 illustrated in Fig. 8 produce a logarithmic 
divergence in the phase measure, whereas the same system for the physical dimension d - 2 has a locally hnite phase 
measure near the critical points. Likewise, the singularities will generally be less important for A-wave resonances 
with N large, since what matters is the size of the phase-space dimension D - (N - 2)d. The cautionary remark to 
these general rules of thumb is that degeneracy degree d > 0 can lead to a stronger singularity at the critical set, and 
this may result in a divergence even when D is larger than 2. This is what occurs in the cases of three-dimensional 
inertial-waves and electron-hole excitations in two-dimensional graphene, for example. 

We collect our case studies in the table below. The table shows space dimension d, order of resonance N, geometry 
of the critical set, degree of degeneracy 6, effective dimensionality D' = d(N- 1)-5 of the phase space, local finiteness 
of the phase measure at the singularity (if any), and divergence or not of the standard collision integral. Note that we 
consider only genuine critical points, not pseudo-critical points (which are usually harmless). When the critical set is 
empty, we take 5 = 0 in the definition of D'. If the phase measure is locally infinite near the singularity, we indicate 
the nature of the divergence: 
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Wave system 

d 

N 

critical set 

6 

D' 

phase measure 

collision integral 

Isotropic power-law, a > 1 

> 2 

3 

0 

- 

d 

finite 

finite 

Acoustic waves 

3 

3 

line 

1 

2 

ill-defined 

divergent 

Rossby/drift waves 

2 

3 

point 

0 

2 

log-divergent 

finite 

Inertial waves 

3 

3 

line 

1 

2 

log-divergent 

divergent 

Internal gravity waves 

3 

3 

0 

- 

3 

finite 

finite 

Surface gravity waves 

2 

4 

0 

- 

4 

finite 

finite 

Surface gravity-capillary waves 

2 

4 

point 

0 

4 

finite 

finite 

Light waves in optical fiber 

1 

4 

point 

1 

1 

linear divergent 

finite (but large) 

Electrons & holes in graphene 

2 

4 

surface 

2 

2 

log-divergent 

divergent 


The diligent reader will recognize that the table is a simplification of the discussion in the text and glosses over some 
of the finer points. (For example, a degenerate critical point with <5 = 1 is possible for gravity-capillary waves at the 
inflection point of the dispersion relation when k — k^.) Nevertheless, the results presented in the table support the 
general lessons educed above. In particular, there tend to be serious consequences for the standard kinetic description 
when the critical set is non-empty and D' < 2. In such cases, closer examination of the problem is warranted, to 
see whether there are any ameliorating circumstances (vanishing interaction coefficients, cancellations in the collision 
integral, etc.) or whether the standard kinetic equation indeed breaks down. Resonance Van Hove singularities and 
their potential effects should be generally recognized as a possibility in wave kinetics. 
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Appendix A. Perturbative Derivation of the Quantum Boltzmann Equation 


The Heisenberg equation of motion for the Hamiltonian (311, ([32) are 


= Ws(^)r™(k) + - 


YjTj f + k2 - ka - k4) 

b i2J3.54^ 

><[7'k:rk"krL,(k2)r.3^(k3)r.4«(k4) - 
We shall write this using a self-explanatory shorthand notation as 

ijla = WlTla + I X J"^^234 <^34 [7’432irlir3/7r4a “ 7’4312rLT3ar4&] ■ 

Hermiticity of the interaction Hamiltonian requires 

r * ^ 

1234 “ f 4321 

and, using the canonical anti-commutation relations, one can also impose the symmetry 


(A.l) 


(A.2) 


Ti234 - T2m- 


One can check that these relations are satisfied for the coefficient T 1234 arising from the Coulomb interaction between 
electrons & holes in graphene by means of the explicit expression given in ca, Eq.(3.9). Using these symmetries, 
we have 

iyia = WlTlfl ^234 ^34 [ 7 ’i 234 rLT 3 fer 4 fl + (3 <-^ 4)] . (A.3) 

Introducing yla = e'^'^yi, we have 


ijjla = I X / *^234 5346 '"^^ [T\234ylhy2b74a + (3 ^ 4)] . 


(A.4) 


Here W 34 = Usiiki) + a)s^{k 2 ) - cjs^ik^) - a>s^(k 4 ). For simplicity, we shall omit the tilde “ ” from now on. We shall 
calculate the mean occupation number rii perturbatively by expanding the creation/annihilation operators into a power 
series 

y\a = Tia + + a^fa + oia^)- (A.5) 

A straightforward calculation gives 


(0) 

yj = 

yia(0) 

(A.6) 

y(') = 

/ Ifl 

- 2 XI J" ^34[7’l234T21,T3fer4a + (3 4)j . 

(A.7) 

4 la 

4 Xi r '^234567 ^34^1234 

^ b,c ^ 


X 

EiiCL>245'’ ^34^^ei'^256iy 2 i,y(,i.y^‘^y2’by4a ~ ^/(^467’ ^ 34 )^35(ny 2 byScyb^yTby4a 


- 

< 4734 ) 5^7 7’;5g7y|^73/,yJ^y6cr7a + (6 7) H- (3 <-4 4) -(- (3 4,6 7) 

(A.8) 


where we employ the standard definitions Il29l : 

A,(x) - I exp(ixs)ds, Et(x,y) - I Aj(x - y) exp(/sy)£/s. (A.9) 

Jo Jo 
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We shall omit (O)-superscript below, when there is no possibility of confusion. The mean occupation number is 
obtained perturbatively by substituting (|a.5|i into 


(rLria) = {yfaJu) + “((rufr'/j) + c.c) + a" ((rr^r®) + c.c + (rufr'/j)) ■ 


(A. 10) 


We assume as initial condition a fermionic quasi-free state with the 2nd-order correlations (yp^yifc) = 
(yvajl/,) andTji = 1 - ni. Substituting (|A.6 i-( A.8 i into (A. 101, we have 




2 /"^234 A,(CU^^) [r*234 + 0 « 4)] 

4 ^f '^234567 ^34^1234 


} 


I 

II 

III 


^[^1567 ( 


E,iajl%- W^2)d^^7’2567 (rLr!irLT5er3ir4fl) 

~ ^^34)<5677’3567 (7'l'at72iT5c'y6c'y7A74a 

- Et{u^^\ w34)^67T4gg7 {^\,^y\i,y3by\yY(>cyii^ 

+ (6 ^ 7) + (3 4) + (3 ^ 4,6 ^ 7) 

^ J" ^^234567 534^567^1'234Af(-W342)A,(W34) 

) + (6 ^ 7) + (3 4) + (3 ^ 4,6 ^ 7) 


1 

4 


ylyly2byly6cy7a 


Using Wick’s rule, we have 




itdi^Qini 

5j,(I + II + III) 


where 


I 

II 

III 

IV 


^ / 


J" £^234 ^34 £’/(0;W34)f?i234 nin3n4+ £',(0;0) J" d2 ^liNT 1221 - Ti2\2)n2 

J" c /234 ^34 Et(0; cl>j^)Ri234 (nin2«4 + nin2«3) “ E,(0; 0) J" c /2 
J" c/234 ^34 Et(0; cl>j^)Ri234 (nin2«3 + nin2«4) - E,(0; 0) J" 


Q [ + ^^2 


(A. 11) 


(A. 12) 

(A. 13) 

(A. 14) 
(A. 15) 
(A. 16) 


(A7’i 221 - 7’i2i2)n2 


2 

Qi + ^2 

c /2 -Z-(A7’i 221 - 7’i2i2)n2 


^234 ^34 |A,(( 7734)1 /?1234 n2«3«4 + ^lA/CO)! 


(A. 17) 


with 


/^1234 - -jTl234{NT\234 - T 1243 ) + (3 4) 

Now, recalling the standard relations ||29]| 

f2 

£,(0;0)=-, A,(0) = f, 


(A. 18) 
(A. 19) 


(A.20) 
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we see that the terms which involve these factors in the previous expressions are undesirable. Their growth is 0{t^) 
at long times, by far the most secular behavior, but they do not correspond to terms in the expected kinetic equation. 
Fortunately, it is straightforward to check that the sum of all these undesirable terms arising from I-IV exactly cancel. 
Also, the 0(a) term from (A. 141 has vanishing real part and thus gives a zero contribution to the evolution of the 
occupation numbers. Then, using 


and 


2Re£',(0;.r) = |A,(x)P = —— 


n\nT,n\ — nin2«4 ~ nin2«3 + niniti'i + nin2«4 ~ n2«3«4 = nin2«3«4 ~ 11111211^114, 


together with (y|,^yia) = we can combine (A. 15 i-( A.16 1 to get 


2 r ^2(1 - cos(ai\^t)) __ 

ni(f) = ni(0) - a t/234 ^34- 7-p^ -R1234 (nin2«3n4 - «i«2«3n4) ■ 

J ("34)^ 


(A.21) 

(A.22) 

(A.23) 


This is equivalent to the expression (33 1 in the text. 

There is another approach in the literature for dealing with the undesirable 0{fi) terms which we should briefly 
mention. It is possible to exactly remove those terms by frequency renormalization umiii, introducing 

u)\ — CJ\ + a^i 


into the Heisenberg equations (A.3 1 , which then becomes 

;yia = OiJia + I ?/ ^234 ^34 [^’1234 {jlb73b74a - {7lb73h)74a + {7lb74a)73b) + (3 « 4)] . (A.24) 

The rest of the derivation is as before, except that now one defines yia = e“^'‘yi using the renormalized frequency. 
The counterterms which appear in the renormalized Heisenberg equations of motion can be readily checked to cancel 
all of the 0(t^) terms in each of the individual expressions I-IV, without the necessity of adding them together. The 


final result is the same as (A.231 except that is replaced with W 34 - cb i + 0)2 - d >2 - 0 ) 4 . 


This alternate procedure yields the same kinetic equation as discussed in the text, but with bare frequencies re¬ 
placed by renormalized frequencies. There might thus naively appear to be an inconsistency between the two ap¬ 
proaches. In particular, the quantum Boltzmann equation for electrons in graphene obtained by the alternate procedure 
would disagree with that derived previously IfTrUTSll . with the bare frequency Wj(k) = svfk undergoing an additional 
renormalization or, equivalently, with an additional renormalization of the Fermi velocity vf- However, this incon¬ 
sistency is only apparent. As is well-known in the wave turbulence literature (e.g. 121 p.71), the two approaches 
lead to equivalent kinetic equations to order 0{a^), since the frequency renormalization is 0{a) and thus corrects the 
kinetic equation only to order 0{a^). The frequency renormalization is therefore entirely optional in the derivation of 
the kinetic equation at order O(a^). The consistency of the two approaches is further evidenced by the fact that the 
undesirable O(t^) terms cancel completely at order O(a^) without any use of a frequency renormalization. 
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